Search for: All records

Novel traits in the order Lepidoptera include prolegs in the abdomen of larvae, scales, and eyespot and band color patterns in the wings of adults. We review recent work that investigates the developmental origin and diversification of these four traits from a gene-regulatory network (GRN) perspective. While prolegs and eyespots appear to derive from distinct ancestral GRNs co-opted to novel body regions, scales derive from in situ modifications of a sensory bristle GRN. The origin of the basal and central symmetry systems of bands on the wing is associated with the expression of theWntAgene in those regions, whereas the more marginal bands depend on two other genes,Distal-lessandspalt. Finally, several genes have been discovered that play important roles in regulating background wing color, via the regulation of pigmentation GRNs. The identification of shared and novelcis-regulatory elements of genes belonging to these distinct GRNs helps trace the developmental and evolutionary history of these traits. Future work should examine the extent to which ancestral GRNs are co-opted/modified to produce the novel traits and how these GRNs map to specific cell types in ancestral and derived traits. 

Evolutionary variation in the wing pigmentation of butterflies and moths offers striking examples of adaptation by crypsis and mimicry. Thecortexlocus has been independently mapped as the locus controlling color polymorphisms in 15 lepidopteran species, suggesting that it acts as a genomic hotspot for the diversification of wing patterns, but functional validation through protein-coding knockouts has proven difficult to obtain. Our study unveils the role of a long noncoding RNA (lncRNA) which we nameivory, transcribed from thecortexlocus, in modulating color patterning in butterflies. Strikingly,ivoryexpression prefigures most melanic patterns during pupal development, suggesting an early developmental role in specifying scale identity. To test this, we generated CRISPR mosaic knock-outs in five nymphalid butterfly species and show thativorymutagenesis yields transformations of dark pigmented scales into white or light-colored scales. Genotyping ofVanessa carduigermline mutants associates these phenotypes to small on-target deletions at the conserved first exon ofivory. In contrast,cortexgermline mutant butterflies with confirmed null alleles lack any wing phenotype and exclude a color patterning role for this adjacent gene. Overall, these results show that a lncRNA gene acts as a master switch of color pattern specification and played key roles in the adaptive diversification of wing patterns in butterflies. 

Abstract Heliconius butterflies, a speciose genus of Müllerian mimics, represent a classic example of an adaptive radiation that includes a range of derived dietary, life history, physiological and neural traits. However, key lineages within the genus, and across the broader Heliconiini tribe, lack genomic resources, limiting our understanding of how adaptive and neutral processes shaped genome evolution during their radiation. Here, we generate highly contiguous genome assemblies for nine Heliconiini, 29 additional reference-assembled genomes, and improve 10 existing assemblies. Altogether, we provide a dataset of annotated genomes for a total of 63 species, including 58 species within the Heliconiini tribe. We use this extensive dataset to generate a robust and dated heliconiine phylogeny, describe major patterns of introgression, explore the evolution of genome architecture, and the genomic basis of key innovations in this enigmatic group, including an assessment of the evolution of putative regulatory regions at the Heliconius stem. Our work illustrates how the increased resolution provided by such dense genomic sampling improves our power to generate and test gene-phenotype hypotheses, and precisely characterize how genomes evolve. 

Abstract Neotropical Heliconius butterflies are well known for their intricate behaviors and multiple instances of incipient speciation. Chemosensing plays a fundamental role in the life history of these groups of butterflies and in the establishment of reproductive isolation. However, chemical communication involves synergistic sensory and accessory functions, and it remains challenging to investigate the molecular mechanisms underlying behavioral differences. Here, we examine the gene expression profiles and genomic divergence of three sensory tissues (antennae, legs, and mouthparts) between sexes (females and males) and life stages (different adult stages) in two hybridizing butterflies, Heliconius melpomene and Heliconius cydno. By integrating comparative transcriptomic and population genomic approaches, we found evidence of widespread gene expression divergence, supporting a crucial role of sensory tissues in the establishment of species barriers. We also show that sensory diversification increases in a manner consistent with evolutionary divergence based on comparison with the more distantly related species Heliconius charithonia. The findings of our study strongly support the unique chemosensory function of antennae in all three species, the importance of the Z chromosome in interspecific divergence, and the nonnegligible role of nonchemosensory genes in the divergence of chemosensory tissues. Collectively, our results provide a genome-wide illustration of diversification in the chemosensory system under incomplete reproductive isolation, revealing strong molecular separation in the early stage of speciation. Here, we provide a unique perspective and relevant view of the genetic architecture (sensory and accessory functions) of chemosensing beyond the classic chemosensory gene families, leading to a better understanding of the magnitude and complexity of molecular changes in sensory tissues that contribute to the establishment of reproductive isolation and speciation. 

Although nearly phenotypically identical, few regulatory regions are conserved between two pairs of butterfly species. 

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis -regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila . We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution. 

Butterfly color patterns provide visible and biodiverse phenotypic readouts of the patterning processes. While the secreted ligand WntA was shown to instruct the color pattern formation in butterflies, its mode of reception remains elusive. Butterfly genomes encode four homologues of the Frizzled-family of Wnt receptors. Here we show that CRISPR mosaic knock-outs of frizzled2 (fz2) phenocopy the color pattern effects of WntA loss-of-function in multiple nymphalids. While WntA mosaic clones result in intermediate patterns of reduced size, fz2 clones are cell-autonomous, consistent with a morphogen function. Shifts in expression of WntA and fz2 in WntA crispant pupae show that they are under positive and negative feedback, respectively. Fz1 is required for Wnt-independent planar cell polarity (PCP) in the wing epithelium. Fz3 and Fz4 show phenotypes consistent with Wnt competitive-antagonist functions in vein formation (Fz3 and Fz4), wing margin specification (Fz3), and color patterning in the Discalis and Marginal Band Systems (Fz4). Overall, these data show that the WntA/Frizzled2 morphogen-receptor pair forms a signaling axis that instructs butterfly color patterning, and shed light on the functional diversity of insect Frizzled receptors. 

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual H. charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for several transcription factor binding site candidates with known function in neuron development in Drosophila. We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution. 

Abstract Sexually dimorphic development is responsible for some of the most remarkable phenotypic variation found in nature. Alternative splicing of the transcription factor gene doublesex (dsx) is a highly conserved developmental switch controlling the expression of sex-specific pathways. Here, we leverage sex-specific differences in butterfly wing color pattern to characterize the genetic basis of sexually dimorphic development. We use RNA-seq, immunolocalization, and motif binding site analysis to test specific predictions about the role of dsx in the development of structurally based ultraviolet (UV) wing patterns in Zerene cesonia (Southern Dogface). Unexpectedly, we discover a novel duplication of dsx that shows a sex-specific burst of expression associated with the sexually dimorphic UV coloration. The derived copy consists of a single exon that encodes a DNA binding but no protein-binding domain and has experienced rapid amino-acid divergence. We propose the novel dsx paralog may suppress UV scale differentiation in females, which is supported by an excess of Dsx-binding sites at cytoskeletal and chitin-related genes with sex-biased expression. These findings illustrate the molecular flexibility of the dsx gene in mediating the differentiation of secondary sexual characteristics.